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Summary 

The analysis procedure of an existing computer pro- 
gram for predicting the off-design performance of radial- 
inflow turbines was extensively modified to improve its 
accuracy and to extend its applicability. The rotor inlet 
slip factor correlation was replaced by one that includes 
rotor blade sweep. At the rotor exit a sector analysis was 
used in place of the original mean-line analysis to model 
the large radial variations occurring there. Disk friction, 
clearance, and vaneless space losses were added to the 
loss model, and the method of computing trailing-edge 
loss was changed. Finally, the entire analysis procedure 
was recoded. 

The resulting computer program, called RTOD (Radial 
Turbine Off-Design), predicts the performance of a 
radial turbine (with or without rotor blade sweep) as a 
function of pressure ratio, speed, and stator setting. The 
program models the flow with stator viscous and trailing- 
edge losses; a vaneless space loss between the stator and 
the rotor; and rotor incidence, viscous, trailing-edge, 
clearance, and disk friction losses. The stator and rotor 
viscous losses each represent the combined effects of 
profile, endwall, and secondary flow losses. The stator 
inlet and exit and rotor inlet flows are modeled by a 
mean-line analysis, but a sector analysis is used at the 
rotor exit. For a pivoting stator the leakage flow through 
the clearance gap is accounted for . 

This report describes the program modeling changes 
and, for completeness, presents all analysis equations. 
The program input and output, including error messages, 
are also described and are illustrated with an example 
case. 

Comparing calculated results with experimental data 
shows that the revised program predicts mass flow rate 
better than does the original program. Calculated results 
are also used to illustrate the potential improvement in 
off-design performance offered by rotor back-sweep for 
high-work-factor radial turbines. 


Introduction 

The radial turbine has been used extensively in 
automotive turbochargers and aircraft auxiliary power 
units and has recently been given serious consideration 
for primary powerplant application for both automobiles 
and rotorcraft. In automotive and rotorcraft applications 
the turbine operates over a wide range of power settings 
and requires accurate off-design performance prediction 
for a successful design. 

A computer code for the off-design performance anal- 
ysis of radial-inflow turbines is described in reference 1 , 
and additions to this program for variable-area stators 
are reported in reference 2. Although this program has 


been used extensively, it still requires modifications to 
extend its accuracy and applicability. Current interest in 
highly loaded radial turbines has resulted in designs 
having rotor back-sweep, a feature that this program can- 
not accommodate. Also, the analysis model, as shown in 
reference 1, predicts too low a value for mass flow in the 
region of rotor choke. In addition, the loss model does 
not include all of the losses present in the design program 
of reference 3, and the trailing-edge loss does not 
properly account for the effect of variable stator angle. 

To overcome these deficiencies, extensive modifica- 
tions were made to the analysis model. The rotor slip 
factor equation was replaced by one that includes rotor 
blade sweep. A sector analysis, similar to that of 
reference 3, was used at the rotor exit in place of the 
mean-line calculation. Disk friction, clearance, and 
vaneless space losses were added to the loss model, and 
the stator and rotor trailing-edge loss equations were 
changed. Finally, the entire procedure was recoded. 

The resulting computer program, called RTOD (Radial 
Turbine Off-Design), which is described in this report, 
computes off-design performance of a radial turbine by 
modeling the flow with stator viscous and trailing-edge 
losses; with a vaneless space loss between the stator and 
the rotor; and with rotor incidence, viscous, clearance, 
trailing-edge, and disk friction losses. RTOD can analyze 
a variable stator, including the leakage flow in the 
clearance gap. Furthermore, RTOD determines the rotor 
exit flow variation from hub to tip and can account for 
the effect of swept rotor blading. The user must specify 
the stator and rotor viscous losses (which represent the 
combined effects of profile, endwall, and secondary flow 
loss) to match a reference performance point. RTOD 
then calculates the turbine performance over a specified 
operating range. 

For the sake of completeness this report describes all of 
the modeling assumptions, but with emphasis on recent 
modifications, comprised in the computer program 
RTOD. The program input and output are described and 
the use of the program is illustrated by example 
problems. Results obtained from this program and from 
the program of reference 1 are compared with experi- 
mental data. RTOD can be obtained from the Computer 
Software Management Information Center (COSMIC), 
University of Georgia, Athens, Georgia 30602. 


Analysis Modification 

This section describes the changes that were made to 
the previous analysis method and loss model of 
references 1 and 2. The described losses are applied at or 
between the flow stations shown in figure 1 . Station 0 is 
the station upstream of the stator, station 1 is just inside 
the stator trailing edge, and station 2 is immediately 


downstream of the stator trailing edge. Station 3 is 
upstream of the rotor, station 4 is just inside the rotor 
trailing edge, and station 5 is immediately downstream of 
the rotor trailing edge. 

Rotor Analysis 

The rotor analysis was extended to include a slip factor 
correlation that can account for swept rotor blading. 
Also, at the rotor exit a sector analysis replaces the 
previous mean -line analysis. The details are as follows 
(all symbols defined in appendix A): 

Slip factor. - In reference 1 the optimum rotor inlet 
tangential velocity is calculated from the Stanitz slip 
factor correlation 


obtained from the velocity diagram work by subtracting 
the losses due to disk friction and tip clearance. This 
means that the sector velocity diagrams are not consistent 
with the shaft output work. 

Loss Model 

In reference 1 the losses through the stator and the 
rotor are each modeled by a single overall loss. In RTOD 
a vaneless space loss (between stations 2 and 3) was 
added, as well as a rotor disk friction and tip clearance 
loss. Furthermore, the method of calculating trailing- 
edge pressure drop was changed. These changes are 
discussed in this section. 

Vaneless space loss. - The program of reference 1 
assumes no pressure drop from the stator exit to the rotor 
inlet (between stations 2 and 3). In RTOD the total 
pressure drop is obtained from 


The Stanitz correlation is limited to radial blading only. 
RTOD uses the Wiesner slip factor correlation of 
reference 4, which is a function of blading (sweep) angle. 
The correlation, when adapted to a turbine geometry, 
takes the form 


Ap' = 


4JlpV2 

2D^c 


where / is the Fanning friction factor for parallel walls 
and / is a length based on an average flow angle between 


^w,3,opt ~ ^3 




-Vcos «bl/(^bl)° '^][l -[(^4,ay/^3-f|im)Al -f|im)J^] 
1 - tan a ^[/ tan 03 


) 


For radial-bladed rotors (no blade sweep) the difference 
between the two correlations is quite small. For the 
example 1 turbine the Wiesner correlation predicts an 
optimum inlet tangential velocity 2.7 percent less than 
that predicted by the Stanitz correlation. The full 
incidence loss calculation procedure is shown by 
equations (B71) to (B82). 

Exit sectors. - The program of reference 1 analyses the 
rotor exit at a mean radius only. In RTOD the rotor exit 
is divided into a specified number of sectors that are 
coupled by simple radial equilibrium 

dp^pV} 

dr rg^ 

Temperature, pressure, and velocity are calculated for 
each sector (eqs. (B94) to (B106)), and overall values of 
mass flow and work are obtained by summing the indivi- 
dual sector values. These calculations take into account 
the viscous loss, since a viscous rotor loss coefficient is 
specified for each sector. The disk friction and tip 
clearance losses, however, are of an overall nature and do 
not reflect hub-to-tip variations. The shaft output work is 


stator outlet and rotor inlet. The calculation procedure is 
given by equations (B45) to (B54). 

Rotor disk friction loss. - The program of reference 1 
does not calculate a separate disk friction loss but lumps 
all losses (viscous, disk friction, clearance, and trailing 
edge) into one overall loss coefficient. In RTOD the 
following expression from reference 3 is used for the disk 
friction loss (eq. (B129)): 

0.02125 pjUj/j 

Rotor tip clearance loss, - The following expression 
from reference 3 is used in RTOD to calculate the rotor 
tip clearance loss (eq. (B130): 

2A^j_5(/^d/2^4,tip) 

“ ^4,hub^^4,tip) 

The formulation assumes that the fractional loss due to 
clearance is equal to the ratio of clearance to passage 
height at the rotor exit. 
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station 



Figure 1. - Radial turbine schematic and station identi- 
fication. 


Trailing-edge pressure drop* - Reference 1 assumes no 
loss in total pressure across the stator or rotor trailing 
edges (stations 1 to 2 and 4 to 5, respectively). This is 
contrary to observed results, since it is known that a total 
pressure drop does occur, although it is usually small. 
Reference 2 assumes that static pressure remains constant 
across trailing edges. Although this assumption gives a 
drop in total pressure, the predicted total pressure drop 
decreases as the flow angle is increased. This is again 
contrary to observed results and will cause errors when 
analyzing a pivoting stator. The expressions for total 
pressure drop across the stator and rotor trailing edges 
used in RTOD were obtained from the expression for 
one-dimensional, compressible fluid flow with sudden 
area increase. The assumed expressions, although not 
based on rigorous analysis, predict the correct trend of 
total pressure loss with flow angle. Furthermore, the 
predicted levels of loss appear reasonable. For the stator 
the drop in total pressure is expressed by (eq. (B39)) 



and for the rotor it is (eq. (B1 13)) 


Program Input 

The program input consists of a title card, the gas 
properties, the turbine geometry description, and the 
parameters that specify the desired calculation options. 
Except for the title card, all data input is in NAMELIST 
format. Data may be input in either SI or U.S. customary 
units. 

Figure 2 identifies the angle and the angle sign conven- 
tion used in RTOD. At the stator and the rotor inlet, all 
angles are defined from a radial line; at the rotor exit, 
angles are defined from an axial line. Figure 3 is a graphic 
representation of the geometry input variables, and 
figure 4 shows the example 1 data input. 

All input variables are listed in table I, along with the 
required units, the default values, and the type of 
variables needed (real or integer). 

The user must specify the viscous stator and viscous 
rotor losses (PRIO or EBR3D for the stator and XKR for 





Figure 2. - RTOD angle sign convention. 
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the rotor) such that the RTOD prediction matches two 
performance parameters (usually flow and work or flow 
and efficiency) at a particular pressure ratio. The 
parameters to be matched can either be design-point 
values obtained from a design-point program or test 
data. 


Program Output 

The program output is in either SI or U.S. customary 
units, corresponding to the chosen input units. The 
output consists of a printout of the title card, all input 
data, and one of three chosen levels of output as 
described in this section. The RTOD data printout 
corresponding to the input data of figure 4 is shown in 
figure 5. 

Short Output 

The short output consists of a tabulation of 12 overall 
parameters for each calculated operating point. Figure 6 
shows the short output corresponding to the input data of 
figure 4. Note that the rotor chokes and that the value of 
stator exit critical velocity ratio (V/ at rotor choke is 


printed out. Also, the rotor limit load point is identified. 
The parameters in the short output are as follows (with 
the corresponding equations in parentheses): 

(1) Equivalent flow (B148) 

(2) Equivalent work (B150) 

(3) Actual flow (B44) 

(4) Actual work (B131) 

(5) Equivalent power (B151) 

(6) Equivalent torque (B 152) 

(7) Equivalent /?o //?5 (B153) 

(8) Equivalent /?o //?5 (B 154) 

(9) Blade-jet speed ratio (B162) 

(10) Work factor (B161) 

(11) Static efficiency (B 134) 

(12) Total efficiency (B135) 

Medium Output 

The medium-output option gives the same tabulation 
of overall parameters as the short output and also gives 
tabulated velocity, temperature, pressure, and flow angle 
information at stations 0 to 5. At stations 3 to 5 the 
information is presented for both the absolute and 
relative frames of reference, and at stations 4 and 5 (rotor 
exit) the information is presented for the hub, mid, and 
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tip sectors only. If an even number of sectors is specified, 
the midpoint sector is taken to be the sector above the 
mean radius. The medium-output option also presents 
additional information for stations 1, 3, and 5. For the 
most part this information is self-explanatory, except for 
these items: 

Station J 


stream flow (pV)\ (equal to 1.0 if 
p{/po was specified; less than 1.0 if 
^30 was specified). 

BLl blockage, defined by ratio of vane 

trailing-edge tangential thickness to 
vane trailing-edge spacing 

Station 3 


EBR3D 


PIT/POT 


V/VCR AT 
MAX FLOW 


three-dimensional kinetic energy loss 
coefficient, ^30, as specified; or, if 
total pressure ratio was specified, 
calculated from equation (B3) 
total pressure ratio, p{/ P q (input 
parameter PR 10), as specified; or if ^30 
was specified, calculated from equation 
(B4) 

critical velocity ratio (F/L^.j-)| that 
produces maximum value of main- 


REL, T. PR. relative total pressure after incidence 
AFTER INCID. loss, a fictitious pressure that shows the 
LOSS effect of incidence 

U3 rotor inlet tip speed 


Figure 7 shows the medium output for a of 0.6 

only (for the example input of fig. 4, except that 
OPT = 1). The output shown is produced for every value 
of {y/Vcr)^ its specified range. 


TABLE I. - INPUT VARIABLES 


Input 

parameter 

Description and units 

Variable 

type^ 

Oefau 1 1 
val ue 

Overall input variables 

TITLE 

Title card (up to 80 alphanumeric characters) 

- 

- 


lUNTS 

Type of units (0 for U.S. customary, 1 for SI) 



0 

TOT 

Total temperature at station 0, K (®R) 

R 


POT 

Total pressure at station 0, N/cm^ (psia) 




G 

Specific heat ratio 




R 

Gas constant, J/kg-K {ft-lbf/lbm-®R) 




XMU 

Gas viscosity, kg/m-s (Ibm/ft-s) 




OPT 

Output option: 

0 - short output 

1 - medium output 

2 - long output 



0 

Stator input variables 

NSTV 

Number of stator vanes 




CLFRl 

Stator total clearance (hub and tip) as fraction of passage height at 
trailing edge 



0.0 

ALO 

Flow inlet angle at station 0, deg 




ALl 

Stator exit blade angle at station 1, deg 




RO 

Radius at station 0, cm (in.) 




R1 

Radius at station 1, cm (in.) 




HSO 

Passage height (including clearance) at stator inlet, cm (in.) 



— 

HSl 

Passage height (including clearance) at stator exit, cm (in.) 




THl 

Vane trailing-edge normal thickness, cm (in.) 




lOPTS 

Option for stator loss calculation: 


I 

1 


1 - user specifies P]/Pq 

2 - user specifies 'e 3 Q 




PRIO 

Specified ratio of P]/Pq (input only if lOPTS = 1) 


R 

0.98 

E8R30 

Specified value of 630 (input only if lOPTS = 2) 



0.055 

VRSTRT 

Initial value of (V/V ) at station 1 (must be subsonic) 
1 




VR5T0P 

Final value of (V/V ) (must be greater than VRSTRT) 
1 




OELVR 

Incremental change in (V/V ) between VRSTRT and VRSTOP 
1 

1 
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TABLE I. - Concluded. 


Input 

parameter 

Description and units 

Variable 

type^ 

Default 

value 

Rotor input variables 

NR8L 

Number of full rotor blades extending from inlet to output 


R 


ISPL 

Splitter blades (0 for no splitters, 1 for splitters) 


I 

0 

RPM 

Rotative speed, rpm 




ALR3 

Rotor inducer sweep angle (0.0 for radial blade), deg 



0.0 

R3 

Radius at station 3, cm (in.) 




R8W0 

Rotor "b“ width at station 3, cm (in.) 




R4HUB 

Rotor exit hub radius, cm (in.) 




R4TIP 

Rotor exit tip radius, cm (in.) 




CL 4 

Clearance height at rotor exit, cm (in.) 




NSCT4 

Specified number of equal height sectors at rotor exit (25 max) 




0PBTA4 

Rotor exit blade angle distribution option 

0 - user specifies exit blade angle for each sector 

1 - program linearly interpolates for exit blade angle at each sector 

from input BTAHUB and BTATIP 



1 

8ETA4 

Rotor exit blade angle for each sector (hub to tip; input only if 
0P8TA4 = 0). deg 

R 


BTAHU8 

Rotor exit blade angle at hub (input only if 0PBTA4 = 1), deg 

f 



BTATIP 

Rotor exit blade angle at tip (input only if 0P8TA4 = 1), deg 

P 



0PTH4 

Rotor exit normal thickness distribution option; 

0 - user specifies normal thickness for each sector 

1 - program linearly interpolates for normal thickness at each sector 

from input TH4HU8 and TH4TIP 



1 

TH4 

Rotor exit normal thickness for each sector (hub to tip; input 
only if 0PTH4 = 0), cm (in.) 

R 



TH4HUB 

Rotor exit normal thickness at hub (input only if 0PTH4 = 1), cm (In.) 




TH4TIP 

Rotor exit normal thickness at tip (input only if 0PTH4 = 1), cm (in.) 




XKR 

Rotor viscous loss coefficient for each sector (hub to tip) 




PNEGI 

Power in incidence loss expression (for negative incidence) 



2.50 

PROS] 

Power in incidence loss expression (for positive incidence) 



1 .75 

P4RDPC 

Incremental percentage reduction in pressure at rotor exit after 
rotor choke 



10.0 


denotes integer; R denotes real. 


4.59 INCH RADIAL TURBINE IN AIR 

4HATT T0T=518.7p P0T=14,696» G=1.4r R=53.35» XHU=1.2UE-5» 

NSTV^14.» AL0=45.p AL1=72.0» R0=2.88B» Rl=2,345» HS0=.5771» HSl=.5771f 
I0FTS=1 p PR10=.987f VRSTRT=.30, VRSTOP=1.20» DELyR=,05. TH1=, 02293, 

NR&L=11.t I3PL=1t RPH-29550,r R3=2.295, RBWD=.5771, R4HUB=.572» R4TIP=1,646, 
CL4=.0145f NSCT4=11, 0PBTA4 = 1» BTAHUB = -32.50f BTAT IP = -56 . 70 , 0PTH4=lr 
TH4HUB=.077, TH4T IP=^ . 03B , XKR = ll<,240r 0PT = 0 iEND 


Hgure 4. - Example 1 data input. 


Long Output 

The long output is similar to the medium output, 
except that for stations 4 and 5 the tabulated inputs are 
presented for each sector. Also, additional rotor loss 
information is presented. For the most part this 
information is again self-explanatory, except for the 
following items: 

EQUIVALENT SPEED defined by equation (B149) 
SPECIFIC DIAM. specific diameter; defined by 

equation (B160) 

SPECIFIC SPEED defined by equation (BI59) 


The long output also presents a further tabulation of the 
following items for each sector at stations 4 and 5: 


Station 4 


R4 

BL4 

W/WCR AT 


radius associated with midpoint of 
sector 

blockage for each sector; defined by 
ratio of free flow area at station 4 to 
free flow area at station 5 
value of (I^/B^cr )4 that produces 
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4.59 INCH RADIAL TURBINE IN AIR 


IUNTS= 0 


f 

T0T=0.51B7E 03 

R . 

POT=0.1470E 02 

FSIA 

G=0.1400E 01 


t 

XMU=0. 1216E-04 LBM/FT*S. 

R=0.5335E 02 FT*LBF/ ( LBh*R ) 

OF*T= 0 


t 





NSTV=0.1400E 02 


7 

CLFRi= 0.0000 




AL0=0.4500E 02 


D£G» 

R0=O.2988E 01 

IN. 

HS0=0.5771E 00 

IN. 

AL1=0.7200E 02 


DEG» 

R1-0.2345E 01 

IN. 

H51=0.5771E 00 

IN. 

TH1=0.2293E-01 


INf 





IORTS= 1 


f 

PR10=0.9370E 00 

. 

EBR3D= 0.0000 

» 

VRSTRT=0.3000E 00 


f 

VRST0P=0. 1200E 01 


DELVR=0.5000E“01 


NRBL^O.llOOE 02 


f 

ISPL= 1 

. 

RPM=0,2955E 05 

RPM. 

ALR3= 0.0000 


DEGf 

R3=0.2295E 01 

IN. 

RBWD=0,5771E 00 

IN. 

R4HUB=0.5720E 00 


INt 

R4TIP=0. 1646E 01 

IN. 

CL4=0. 1450E-01 

IN. 

NSCT4= 11 


r 





0FTH4= 1 


> 





TH4HUB=0.7700E-01 


IN* 

TH4TIP=0.3800E-01 

IN. 



0PBTA4= 1 


* 





BTAHUB=-.3250E 02 


DEG. 

BTATIF-- .5670E 02 

DEG. 



PNEGI^0.2500E 01 


f 





PPOSI=0. 1750E 01 







P4RDPC=0. lOOOE 02 







BETA4( 1)=-.3360E 

02 

DEG. 

TH4( 1)=0.7523E-01 

IN. 

XKR( 1)=0.2400E 00 


( 2)=-.3580E 

02 


( 2)=0.7168E-01 


( 2)=0.2400E 00 


( 3)=-.3800E 

02 


( 3)=0.6814E-01 


( 3)=0.2400E 00 


( 4)=-.4020E 

02 


( 4)=0.6459E-01 


( 4)=0.2400E 00 


( 5)==-.4240E 

02 


{ 5)=0.6105E-01 


{ 5)=0.2400E 00 


( 6)=-.4460E 

02 


( 6)=0.5750E-01 


( 6)=0.2400E 00 


( 7>=-.4680E 

02 


( 7)=0.5395E-01 


( 7)=0.2400E 00 


( 8)=-.4900E 

02 


< 8)=0.5041E-01 


( 8)=0.2400E 00 


{ 9)=-.5120E 

02 


( 9)=0.46B6E-Oi 


( 9)=0.2400E 00 


(10)=-.5340E 

02 


(10)=0,4332E-01 


(10)=0.2400E 00 


(11)=-.5560E 

02 


(11)=0.3977E'01 


(11 )=0.2400E 00 



Figure 5. - RTOD data printout {corresponding to fig. 4 data input). 


<V/VCR)1 EQUIV, 

EQUIV. 

ACTUAL 

ACTUAL 

EQUIV 

FLOW 

WORK 

FLOW 

WORK 

POWER 

(LBM/ 

(BTU/ 

(LBM/ 

(BTU/ 

(HP) 

SEC) 

LBM) 

SEC) 

LBM) 



0.3000 

0.3769 

4.6418 

0.3769 

4.6421 

2.4756 

0.3500 

0.4337 

6.1643 

0.4337 

6.1646 

3.7830 

0.4000 

0.4878 

7.6727 

0.4B78 

7.6732 

5.2962 

0.4500 

0,5389 

9.1776 

0.5388 

9.1782 

6,9979 

0.5000 

0.5865 

10.6883 

0.5B65 

10.6890 

8.8703 

0.5500 

0.6306 

12.2155 

0.6305 

12.2162 

10.8994 

0.6000 

0.6707 

13.7671 

0.6706 

13.7679 

13.0643 

0.6500 

0.7066 

15.3593 

0,7066 

15.3602 

15.3556 

0.7000 

0.7381 

17.0213 

0.7381 

17.0223 

17.7751 

0.7500 

0.7651 

18.8252 

0.7651 

18.8263 

20.3776 

0,8000 

0.7873 

21 .1710 

0.7B73 

21.1723 

23.5866 

0.8062 

0.7898 

21.9757 

0.7897 

21.9769 

24.5496 


EQUIV. 

EQUIV. 

EQUIV. 

BLADE-JET 

WORK 

STATIC 

TORQUE 

P0T/P5 

P0T/P5T 

SPEED 

FACTOR 

EFFIC. 

(FT*LBF) 



RATIO 




0.4400 

1.2412 

1.2099 

0.9689 

0.3425 

0.6230 

0.6724 

1.2874 

1.2571 

0.6983 

0.4520 

0.7112 

0.9414 

1.3395 

1.3085 

0.8375 

0,5606 

0.7695 

1.2438 

1.3998 

1.3660 

0.7832 

0.6692 

0.8049 

1.5766 

1.4712 

1.4320 

0.7336 

0.7783 

0.8223 

1.9373 

1 * 5575 

1.5090 

0.6875 

0.8886 

0.8255 

2.3221 

1.6623 

1.5989 

0.6448 

1.0008 

0.8185 

2.7293 

1.7930 

1.7058 

0.6048 

1,1160 

0.8032 

3.1594 

1.9637 

1.6373 

0.5661 

1,2362 

0.7799 

3.6219 

2.2084 

2.0099 

0.5267 

1,3668 

0.7467 

4.1923 

2.7154 

2.3076 

0.4757 

1 .5366 

0.6851 

4.3635 

3.0194 

2.4488 

0.4556 

1.5949 

0.6522 


ROTOR HAS CHOKED AT UlVCRl =0 . 3062E 00 ( 4 ITERATIONS) 

HAy P4(l) FOR CHOKED R0TOR=O . 4 1 82E 01 FSIA ( 14 ITERATIONS) 


0.8062 

0.7898 

22.3802 

0.7897 

22.3814 

25.0014 

0.8062 

0.7898 

22.7041 

0.7897 

22.7055 

25.3634 

0.8062 

0.7898 

22.9340 

0.7897 

22.9353 

25.6201 

0.8062 

0.7898 

22.9930 

0.7897 

22.9943 

25.6860 


ROTOR LIMIT LOADING ACHIEVED AT F'4 ( 1 ) =0 . 29 1 9E 01 


4.4438 

3.2415 

2.5399 

0.4438 

1.6242 

0.6301 

4,5081 

3.4313 

2.6337 

0.4330 

1.6477 

0.6085 

A * 0 u v/ S 

3.7395 

2.7322 

0.4231 

1.6644 

0.5869 

4.5655 

3.8505 

2.7747 

0.4193 

1.6686 

0.5779 


rSIA ( 13 ITERATIONS) 


Figure 6. - RTOD short output for figure 4 data input. 


TOTAL 

EFFIC. 


0.7041 

0.7827 

0.8338 

0*8649 

0.8808 

0.8850 

0.6815 

0.8721 

0.8573 

0,8365 

0.3004 

0.7821 


0.7691 

0.7548 

0.7382 

0.7305 
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C;/UCR) 1 

EQUIV. 

FLOW 

(L&rt/ 

SEC) 

EQUIV. 

WORK 

(EiTU/ 

LBH) 

ACTUAL 

FLOW 

<LRM/ 

SEC) 

ACTUAL 

WORK 

(BTU/ 

LBM) 

EOUIV, 

POWER 

(HP) 

EQUIU. 
TO ROUE 
(FT*LBF) 

EQUIV. 

P0T/P5 

EQUIU. 

P0T/P5T 

BLADE-JET 

SPEED 

RATIO 

WORK 

FACTOR 

STATIC 

EFFIC, 

TOTAL 

EFFIC. 

0.6000 

0.6707 

13.7671 

0.6706 

13.7679 

13.0643 

2.3221 

1.6623 

1 .5989 

0.6448 

1.0008 

0.8185 

0.881! 


STATION 

VCR 

V/UCR 

VU/VCR 

0R(VX)/VCR 

IT 

T 

PT 

F' 

FLOW 



(AND 

OR 

OR 

OR 

DR 





ANGLE 



SECTOR) 

WCR 

W/WCR 

WU/UCR 

UR(WX)/UCR 









(FT/S) 




(R) 

(R) 

(PSIA) 

(PSIA) 

(DEG) 



0 

1019.1802 

0.1693 

0.1197 

0.1197 

518.7000 

516.2205 

14.6960 

14.4516 

45.0000 



1 

1019. 1B02 

0.6000 

0.5706 

0.1854 

518.7000 

487,5779 

14.5050 

11.6806 

71.9999 



') 

1019.1802 

0.5960 

0.5706 

0.1721 

518.7000 

487.9883 

14.4916 

1 1 .7042 

73,2137 



3 

1019. 1B02 

0.6095 

0.5831 

0.1775 

518.7000 

4B6.5S76 

14.4669 

11 .5673 

73.0723 



3(R) 

989.8862 

0. 1327 

0.0025 

0.1827 

489.3105 

486.5876 

11.7954 

11.5673 

0.7704 



4 ( 1 ) 

959.0784 

0.3021 

-0.0326 

0,3003 

459.3274 

452,3423 

9,2163 

8.7350 

-6.1953 



4( 6) 

960.0049 

0.3009 

0.0012 

0.3009 

460.2151 

453,2712 

9.2145 

8.7370 

0,2269 



4(11) 

960.7744 

0.2855 

0.0120 

0.2853 

460.9534 

454.6897 

9.1659 

8.7373 

2.4163 



4 ( 1 ) ( R ) 

962 . 1702 

0.3594 

-0.193? 

0.2993 

462.2937 

452.3423 

9.4262 

8.7350 

-33.6000 



4( 6)(R) 

967,0215 

0.4195 

-0.2945 

0.29B7 

466.9673 

453.2712 

9.6964 

B.7370 

-44.5999 



4 ( 1 1 M R ) 

974.5627 

0.4978 

-0.4108 

0.2312 

474.278S 

454.6897 

10.1273 

8.7373 

^55.5999 



5( 1 ) 

959.0784 

0.2231 

-0.0326 

0.3003 

459,3274 

455.5159 

9.1831 

S.9192 

-8.4009 



5( 6) 

960.0049 

0.2603 

0.0012 

0.3009 

460.2151 

455*0198 

9.2060 

8.8474 

0.2624 



5(11) 

<J60.7744 

0.2622 

0.0120 

0.2853 

460.9534 

455.6709 

9.1631 

8.8008 

2.6314 



STATION 1 












PRIMARY WEIGHT 

FLOW ^ 0.67 

06E 00 < 

(LBM^S) 


SPECIFIED 

PRIMARY 

FLOW ANGLE 

- 0.7200E 

02 

( DEG ) 

CLEARANCE WEIGHT 

FLOW ^ 

0.0000 

(LBM./S) 


ACTUAL 

PRIMARY 

FLOW ANGLE 

- 0.7200E 

02 

(DEG) 

TOTAL WEIGHT 

FLOW : 0.67 

06E 00 ' 

!LBM/S) 


CLEARANCE 

FLOW ANGLE 

- 0,0000 

■DEG) 


EBR3D C . 55 

23E-01 




V/UCR AT MAX FLOW 

^ O.IOOOE 

01 



PIT 

/POT ^ 0.95 

7 0E 00 





ELI 

- 0.7051E- 

01 



STATION 3 

OPTIMUM BETA = -.2060E 02 

(DEG) 

REL. T. PR. AFTER INCID. LOSS 

- 0.1177E 02 

■PSIA 

INCIDENCE ANGLE ^ 0.2137E C2 

(DEG; 

U3 

- 0.591SE 03 

(FT/S 


TAT ION 5 







HASS AVERAGED T5 

^ 0.4553E 

03 MASS AVG. 

DELH (STA. 0 TO 5' 

^ 0.16S2E 

02 

',BTU./L&n 

MASS AVERAGED T5T 

- 0.4604E 

03 MASS AVG. 

DELHI (STA. D TO S) 

- 0. 15c-2E 

02 

BTU/ LBM 

FT, AT THE HUB 

^ 0.E91BE 

0 1 F' 5 1 H ' li M S- S A V G , 

POWER (STA, C TO 5; 

^ 0.14e2E 

02 

(HP > 

P5 AT THE SHROUD 

- O.SSOIE 

Oi (PSIA; MASS AVG, 

TORQUE (STA. 0 TO 5^ 

= 0.26-34E 

01 

(FT#LBF ^ 

AVI RAGE PL 

- 0 .866 0 E 

01 ; F' 5 I A ; 





AVEFhOE ftt 

-- O.'^TllE 

01 







Figure 7. - RTOD medium output at (V/ 

» = 0. 6 for figure 4 data 

input. 




MAX FLOW 
U4 


maximum (choking) value of (pW 04 
(always less than 1 .0) 
rotor sector speed 


The long output shown is produced for every value of 
{V/Vizr)^ in its specified range. 


Station 5 


AEROD. ETA 
AEROD. ETAT 
WM(K)/WMTOT 


sector aerodynamic static efficiency; 
defined by equation (B127) 
sector aerodynamic total efficiency; 
defined by equation (B128) 
fraction of total weight flow passed 
by each sector 


Figure 8 shows the long output for a ( V/ of 0.6 only 
(for the example input of fig. 4, except that OPT = 2). 


Error Messages 

All DO-loops in RTOD have appropriate error 
messages that are printed out if convergence is not 
achieved. The specified convergence criteria are defined 
by the ratio of the absolute value of the new value minus 
the old value to the old value and are usually required to 
be less than 1 x 10“^. These are tight tolerances, but no 
DO-loop has ever failed to converge for physically 
possible flow. 
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(p/ycR)i 

EQUIV. 

FLOW 

(LBM/ 

SEC) 

EQUIV. 

WORK 

(BTU/ 

LBH) 

ACTUAL 

FLOW 

(LBM/ 

SEC) 

ACTUAL 

WORK 

tBTU/ 

LBM) 

EQUIV, 

POWER 

<HP> 

EQUIV. 

TORQUE 

(FT*LBF) 

EQUIV, 

P0T/P5 

EQUIV, 

P0T/P5T 

BLADE-JET 

SPEED 

RATIO 

WORK 

FACTOR 

STATIC 

EFFIC, 

TOTAL 

EFFIC. 

0,6000 

0,6707 

13,7671 

0,6706 

13.7679 

13.0643 

2 , 3221 

1.6623 

1.5989 

0,6448 

1,0008 

0.8185 

0.881! 


STATION 

VCR 

V/VCR 

VU/VCR 

VR( VX)/VCR 

TT 

T 

PT 

P 

PLOW 


(AND 

OR 

OR 

OR 

OR 





ANGLE 


SECTOR) 

WCR 

U/UCR 

WU/UCR 

WR(UX)/WCR 








(FT/5) 




(R) 

(R) 

(PSIA) 

(PSIA) 

(DEG) 


0 

1019.1802 

0.1693 

0.1197 

0.1197 

518.7000 

516.2205 

14,6960 

14.4516 

45.0000 


1 

1019,1802 

0.6000 

0.5706 

0.1854 

518.7000 

487.5779 

14.5050 

1 1 .6806 

71,9998 


2 

1019.1802 

0.5960 

0.5706 

0.1721 

518.7000 

487.9883 

14.4916 

1 1.7042 

73.2136 


3 

1019,1802 

0.6095 

0.5831 

0.1775 

518.7000 

496.5876 

14,4669 

11.5673 

73.0722 


3(R) 

989,8862 

0.1827 

0.0025 

0. 1827 

489.3105 

486.5876 

11.7955 

1 1.5673 

0.7704 


4( 1) 

959.0784 

0.3021 

-0.0326 

0.3003 

459.3274 

452.3423 

9.2163 

8.7350 

-6.1954 


4( 2) 

959,2175 

0.3015 

-0.0237 

0.3006 

459.4604 

452.4976 

9.2159 

8,7362 

-4,5008 


4( 3) 

959 . 3916 

0.3015 

-0.0159 

0.3011 

459.6274 

452.6646 

9.2163 

8,7368 

-3.0167 


4( 4) 

959.5384 

0.3015 

-0,0092 

0.3014 

459,8159 

452.8477 

9.2166 

8.7370 

-1.7400 


4( 5) 

959.7964 

0.3014 

-0.0035 

0,3014 

460.0151 

453.0496 

9.2163 

8.7370 

-0.6631 


4( 6) 

960.0049 

0,3009 

0.0012 

0,3009 

460,2151 

453.2712 

9.2145 

8.7370 

0.2268 


4( 7) 

960.2041 

0.2997 

0.0049 

0,2997 

460.4063 

453.5137 

9.2107 

8.7370 

0.9447 


4( B) 

960.3867 

0.2977 

0.0078 

0,2976 

460.5813 

453.7766 

9.2043 

8.7371 

1 . 5058 


4( 9) 

960.5454 

0.2948 

0.0099 

0,2946 

460.7336 

454.0603 

9.1949 

8.7371 

1 .9278 


4(10) 

960,6760 

0.2908 

0.0113 

0.2906 

460.8509 

454.3647 

9.1822 

8.7372 

2,2257 


4( li) 

960,7744 

0.2855 

0.0120 

0.2853 

460.9534 

454.6897 

9,1659 

8.7373 

2,4162 


4( 1)(R) 

962,1702 

0,3594 

-0.1989 

0.2993 

462.2937 

452.3423 

9.4262 

8.7350 

-33.6000 


4( 2)(R) 

962.9231 

0.3692 

-0.2160 

0.2995 

463.0173 

452.4976 

9,4680 

8.7362 

-35.8000 


4< 3) (R) 

963.7849 

0,3803 

-0.2341 

0.2997 

463.3467 

452.6646 

9,5158 

8.7368 

-38.0000 


4 ( 4 ) ( R ) 

964.7554 

0.3925 

-0.2533 

0.2998 

464,7313 

452.3477 

9,5697 

8.7370 

-40.2000 


4 ( 5 ) ( R ) 

'J65.G345 

0.4056 

-0.2735 

0.2995 

465.8215 

453.0496 

9.6299 

8.7370 

-42 . 4000 


4 ( 6 ) •; R ) 

967,0215 

0.4195 

-0.2945 

0.2937 

4*6.9673 

453.2712 

9.6964 

8.7370 

-44.599? 


4( 7)(R) 

963,3162 

0.4341 

-0.3164 

0.2972 

468.2185 

453.5137 

9.7694 

8.7370 

-46.7999 


4( 8)(R) 

969.7185 

0.4493 

-0.3391 

0.2948 

469.5754 

453.7766 

9.8489 

S .7371 

-48 . 9999 


4( 9)<R) 

971.2268 

0.4650 

-0,3624 

0.2914 

471.0376 

454.0603 

9.9350 

8,7371 

-51 . 1999 


4 ( 1 0 ) ( R ) 

972 .3420 

0 . 4812 

-0,3863 

0.28*9 

472.6055 

454.3647 

10.0277 

8.7372 

-53.3999 


4 n n '! R ) 

974 . 5627 

0.4978 

-0.4108 

0.2812 

474.2708 

454.6897 

10.1273 

8 ,7373 

-55.5998 


5 ( 1) 

959.0784 

0.2231 

-0.0326 

0.3003 

459.3274 

455.5159 

9.1331 

8.9192 

-8 . 4009 


O' 2) 

959.2175 

0.2341 

-0.0237 

0 , 3006 

459,4604 

455.2632 

9 .1920 

8,9015 

-5.8013 



959.39-6 

0.2429 

-0,0159 

0.3011 

459.6274 

455.1062 

9,1984 

8.8856 

-3.7447 


J( 4) 

959,5884 

0.2501 

-0 . 00 9 2 

0.3014 

459.8159 

455.0222 

9.2029 

6,3715 

-2.0980 


5( 5) 

959.79o4 

0.2558 

-0 . 0C35 

0 . 3014 

460.0151 

454.9966 

? ( j 0 5 5 

8. 3588 

-0.7812 


5 ( 6 ) 

9*0.0049 

0.3603 

0.0012 

0,3009 

460.2151 

455.0198 

9.2060 

6.8474 

0.2623 


5( 7) 

■^’6-:).2041 

0,263 4 

0.0049 

0 . 2'=’97 

460.4063 

455. 0840 

9,203? 

8.3368 

1.0751 


8) 

960,3867 

0.2652 

0,0078 

0. 2^^6 

460.5S13 

455. 1S43 

9.1908 

6.8270 

1.6909 


5 (. 9 ) 

960.5454 

0.3656 

0.0099 

0.2946 

460 . 7336 

455.3169 

9.1905 

8.8178 

2.1399 


2(10) 

9 6 0 1 6^60 

0.2646 

0.0113 

0.2906 

4 60. 3 j 3 7 

455.4000 

? . 1766 

8 . 309 1 

2 .4458 


5(11 ) 

960.^^44 

0 .2622 

0.0120 

0 . 2S53 

460.9534 

455.670? 

9.1631 

S. 8008 

2.6313 



Figure 8. - RTOD long output at (V/ 


0.6 for figure 4 data input. 


Main Program 

The user can specify geometries and operating condi- 
tions that have no flow solutions. For low flow rates 
(small specified values of V/ at station 1 ) it may not be 
possible to pass the flow at station 4 while maintaining 
radial equilibrium. In a real turbine the flow would be 
passed but with separated or reversed flow in portions of 
the rotor exit. However, RTOD cannot account for that 
situation, and if it occurs, the following error message is 
printed out: 

STATION 4 HAS NOT CONVERGED IN 
25 ITERATIONS FOR (V/VCR)1 = x.xxxxx ±xx 

ROTOR WEIGHT FLOW = xxxxxx d=xx, 

STATOR WEIGHT FLOW = x.xxxxx ±xx 

After this message, program execution starts over with 


the next specified value of stator exit critical velocity 
ratio. 

Subroutine MAXVAL 

Subroutine MAXVAL calculates the maximum value 
of density times velocity pV and the associated critical 
velocity ratio at the stator and rotor exits (stations I and 
4). The input points must define a curve that shows a 
maximum somewhere in the specified interval (other than 
at the end points). If unrealistically large values of stator 
or rotor loss are specified, the shape of the generated 
curve may not exhibit the required characteristic. If this 
happens, the following error message is printed out: 

WARNING-THE CURVE TO BE MAXIMIZED 
IN SUBROUTINE MAXVAL IS NOT OF THE 
PROPER FORM 
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STATION 1 

PRIHARY WEIGHT FLOW = 0.6706E 00 (LPH/S) 
CLEARANCE WEIGHT FLOW = 0.0000 (LPM/S) 

TOTAL WEIGHT FLOW = 0.6706E 00 (LE<h/5> 

EBR3D = 0.5523E-01 
RIT/POT = 0.9870E 00 


STATION 3 

OPTIHUH BETA = -.2060E 02 (DEG) 
INCIDENCE ANGLE = 0.2137E 02 (DEG) 


STATION 5 

MASS AVERAGED T5 = 0.4553E 03 (R) 

MASS AVERAGED T5T = 0.4604E 03 CR) 

P5 AT THE HUB = 0.8918E 01 (PSIA) 

P5 AT THE SHROUD = 0*8801E 01 (PSIA) 

AVERAGE P5 = 0.8860E 01 (PSIA) 

AVERAGE P5T = 0.9211E 01 (PSIA) 


DELHT ACROSS ROTOR = 0.1400E 02 (BTU/LBM) 
DISK FRICTION LOSS = 0,4306E-01 (BTU/LBM) 
TIP CLEARANCE LOSS ^ 0.1B90E 00 (BTU/LBM) 


EQUIVALENT SPEED = 0.2955E 05 (RPM) 
SPECIFIC SPEED = 0.7071E 00 


SPECIFIED PRIMARY FLOW ANGLE = 0.7200E 02 (DEG) 

ACTUAL PRIMARY FLOW ANGLE = 0.7200E 02 (DEG) 

CLEARANCE FLOW ANGLE = 0,0000 (DEG) 

V/VCR AT MAX FLOW = O.lOOOE 01 

BLl = 0.7051E-01 


REL. T. PR. AFTER INCID, LOSS = 0.1177E 02 (PSIA) 
U3 = 0.5918E 03 (FT/S) 


MASS AVG. DELH (STA. 0 TO 5) = 0.16B2E 02 (BTU/LBM) 
MASS AVG. DELHT (STA. 0 TO 5) = 0.1562E 02 (BTU/LBM) 
MASS AVG. POWER (STA. 0 TO 5) = 0.14B2E 02 (HP) 

MASS AVG. TORDUE(STA. 0 TO 5) = 0.2A34E 01 (FTtLBF) 


ACTUAL DELHT = 0.1377E 02 (BTU/LBM) 

ACTUAL POWER = 0.1306E 02 (HP) 

ACTUAL TORQUE = 0.2322E 01 (FTtLBF) 


SPECIFIC DIAM. = 0.2A77E 01 


STATION 4 STATION 5 



R4 


U/WCR AT 

U4 

AEROD. 

AEROD. 

WM(K)/ 

SECTOR 

(IN) 

BL4 

MAX FLOW 

(FT/S) 

ETA 

ETAT 

WMTOT 


1 

0.6208E 

00 

0.2547E 

00 

0.8653E 

00 

0. UOIE 

03 

0.8608E 

00 

O.9105E 

00 

0.4412E-01 

2 

0.71B5E 

00 

0.2154E 

00 

0.9653E 

00 

0.1853E 

03 

0.8557E 

00 

0.9102E 

00 

0.5381E-01 

3 

0.8U1E 

00 

0, 1855E 

00 

0. 8653E 

00 

0.2104E 

03 

0.8505E 

00 

0.9089E 

00 

0.6354E-01 

4 

0.9137E 

00 

0, 1620E 

00 

0.8653E 

00 

0.2356E 

03 

0.8453E 

00 

0.9069E 

00 

0.7325E-01 

s 

0. lOllE 

01 

0. 1431E 

00 

0.8653E 

00 

0.2608E 

03 

0.B403E 

00 

0.9044E 

00 

0.8289E-01 

6 

0. IK'^E 

01 

0.1275E 

00 

0.8653E 

00 

0.28i0E 

03 

0.8354E 

00 

0.9014E 

00 

0.9237E-01 

7 

0. 1207E 

01 

0. U44E 

00 

0.8655E 

00 

0.3112E 

03 

O.B309E 

00 

0.8990E 

00 

0.1016E 00 

e 

0. 1304E 

01 

0.1031E 

00 

0.8653E 

00 

0.3363E 

03 

0.B267E 

00 

0,8943E 

00 

O.1104E 00 

V 

0, 14«)2E 

01 

0.9340E- 

■01 

0.Si53E 

00 

0.3615E 

03 

0.8230E 

00 

0.B904E 

00 

0.1187E 00 

10 

0. 1500 E 

01 

0.8482E- 

-01 

0.3653E 

00 

0. 3867E 

03 

0.8197E 

00 

0,8862E 

00 

0.1263E 00 

11 

0. 1597E 

01 

0.7717E- 

01 

0,8i53E 

00 

0.4119E 

03 

0.B170E 

00 

0.8B1BE 

00 

0.1331E 00 


Figure 8. - Concluded. 


Program execution continues after this message. The 
value of pV used is the greater of the two at the interval 
end points. The program user should reexamine and 
change the specified stator or rotor losses if this error 
message occurs. 


Example Problems 

A radial turbine that was tested in air at pressure ratios 
higher than choking is analyzed in this section, and the 
results are compared with experimental data, together 
with the results from the program of reference 1. Also, 
the same turbine is analyzed for various values of rotor 
back-sweep to illustrate the potential off-design 
performance improvements in high-work-factor turbines. 


Example 1 - Comparison with Test Data and 
Reference 1 Code Prediction 

Figure 9 shows a comparison of calculated versus 
experimental mass flow at 100- and 60-percent design 
speeds for the 11.66-cm (4.59-in.) diameter turbine of 
reference 5. The test data are previously unpublished air 
data that were used because they covered a wider range of 
pressure ratio (beyond rotor choke) than the published 
argon data. Both the RTOD and reference 1 code predic- 
tions were matched with data at the design equivalent 
total-to-static pressure ratio of 1.54. The RTOD results 
more closely predict the experimental mass flow rates at 
high pressure ratios. At both 100- and 60-percent design 
speeds, RTOD predicts an approximately 1 percent 
greater mass flow ratio than the reference 1 code. The 
improved mass flow prediction is attributed to the more 
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Figure 9. - Mass flow rate comparison. 



accurate geometry and flow modeling at the rotor trailing 
edge. Both the RTOD and reference 1 code total effi- 
ciency predictions agree closely with test data as shown in 
figure 10. 

Although the RTOD mass flow predictions are an im- 
provement over the reference 1 code predictions, they still 
underpredict choking weight flow. At 100-percent design 
speed the RTOD prediction is 1.63 percent below test 
results; while at 60-percent speed it is 2.28 percent low. 
This discrepancy could be caused by the basic rotor loss 
modeling scheme. Rotor viscous loss is taken to be a 
function of rotor relative inlet and exit velocities squared 
(eq. (B87)), multiplied by a loss coefficient (input param- 
eter XKR) that is assumed to remain constant. This as- 
sumption forces the calculated maximum mass flow rate 
per unit area to occur at low values of rotor exit relative 
critical velocity ratio = in example 1). 

Proper variation of the rotor loss coefficient with relative 
flow velocity might improve the choking mass flow 
prediction. 

Example 2 - Effects of Rotor Back-Sweep 

To illustrate the effects of rotor back-sweep, the 
example 1 turbine was analyzed with assumed back- 
sweep angles of \ 5 ° and 30°. The predicted total 
efficiencies are shown in figure 11. For high values of 
work factor (the region where modern radial turbines are 
likely to operate) the falloff in efficiency is less with back- 
sweep than for a radial-bladed rotor (because the 
incidence is lower). 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio, March 25, 1983. 
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Appendix A 
Symbols 


A area, m2 (ft2) 

BL4 ratio of blocked area at station 4 to free flow 

area at station 5 

CLFRI ratio of total clearance area to total passage 
area at station 1 

D/j hydraulic diameter, m (ft); defined by 

eq. (B46) 

specific diameter; defined by eq. (B160) 

^30 three-dimensional kinetic energy loss 

coefficient; defined by eq. (B3) 

/ friction factor 

force-mass conversion constant, 1 (32.174 
(lbm-ft)/(lbf-s2)) 

/i height, m (ft); or enthalpy, J/kg (Btu/Ibm) 

J dimensional constant, 1 (778 ft-lbf/Btu) 

A rotor viscous loss coefficient 

/c sector number 

L kinetic energy loss, J/kg (Btu/lbm) 

/ length, m (ft) 

M Mach number 

TV number of sectors at station 4 

N^\ number of rotor blades (including splitters) 

Ns specific speed; defined by eq. (B159) 

n incidence exponent 

P power, W (hp) 

p pressure, N/m2 (psia) 

Q parameter defined by eq. (B166) 

R gas constant, J/kg-K (ft-lbf/lbm-°R) 

Re Reynolds number; defined by eq. (B47) 

r radius, m (ft) 

T temperature, K (®R) 

TRQ torque, N-m (ft-lbO 

U tangential velocity, m/s (ft/s) 

V velocity in absolute reference frame, m/s 

(ft/s) 

W velocity in relative reference frame, m/s (ft/s) 

WRK work, J/kg (Btu/lbm) 

w flow rate, kg/s (Ibm/s) 

a fluid now angle in absolute frame of refer- 

ence; or blading angle, deg 

0 fluid flow angle in relative frame of reference; 

or blading angle, deg 
7 ratio of specific heats 

6 parameter defined by eq. (B156) 


6 parameter defined by eq. (B157) 

fiim parameter defined by eq. (B73) 

y] efficiency 

^cr parameter defined by eq. (B155) 

/X viscosity, kg/m-s (Ibm/ft-s) 

V loading coefficient (blade-jet speed ratio); 

defined by eq. (B162) 
p density, kg/m^ (Ibm/ft^) 

i' work factor; defined by eq. (B161) 

Q rotational speed, rpm 

Subscripts: 

act actual 

aed aerodynamic 

av average 

bl blading 

ck choked 

cl clearance 

cr critical 

df disk fraction 

equ equivalent 

/ friction 

hub hub 

/ incidence 

id ideal 

mx maximum 

opt optimum 

p primary 

r radial component 

ref reference 

shr shroud 

std standard 

tip tip 

tot total 

u tangential component 

X axial component 

0,1,2, station numbers defined by fig. 1 

3,4,5 

Superscripts; 

( )' total conditions in absolute reference frame 

( )' ' total conditions in relative reference frame 

( ) mass averaged 
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Appendix B 
Equations 


Main Program (MAINP) 

The analysis of the flow through the turbine starts at 
the stator exit (station 1) for each specified value of 
iV! conditions are then calculated at station 0 

and at stations 2 to 5, The procedure is described in this 
appendix. 

Station 1— primary //ow. -Total temperature at 
station 1 is assumed to be equal to that at station 0. Thus 

T[.p=n.p (Bi) 

The maximum value of primary flow per unit area 
[(pK)j and the value of {V/Vqx)x where that 

maximum occurs l(V/Vcr), 1 ^^e calculated by 
subroutine MAXVAL. If the total pressure drop through 
the stator (station 0 to station 1) is modeled by specifying 
the total pressure ratio p{/poy nmx 

(K/Kcr)j ^^=1.0. If the total pressure drop is modeled 
by specifying ^30, it occurs at (l^/l^cr),^nix 
p{/po is specified, total pressure at station 1 is calculated 
from 

Pi,- (B2) 


Both methods of calculating the total pressure drop are 
equivalent to specifying an overall viscous loss that takes 
into account the vane profile, endwall, and secondary 
flow losses. 

From the isentropic relationships 


II 


(B5) 

II 


(B6) 

and 



01 = 


(B7) 

II 


(B8) 


and the three-dimensional kinetic energy loss coefficient For values of {V/ less than (Vi primary 

is obtained from flow exits at the specified stator exit blade angle a\. Thus 



For values of ( V/ Vcr)^ p greater than ( V/ Vqt)^ p (choked 
flow) the exit angle is decreased to maintain maximum 
primary flow rate. The choked-flow exit angle is 


(BIO) 


If e^o is specified, 



The radial and axial components of the primary flow are 
then defined by 

,p “ sin Of] (Bll) 

^r,\,p~ ^l,p ^l,p (B12) 


(B4) 


Mach number is calculated from 
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(B13) 


If the radial component of Mach number (A/j cos aj p) 
is greater than 1 .0, the specified ( V/ is reduced until 
Mj p cos 1.0 (defined as limit loading). 

Station I -clearance flow.-J\\c clearance flow is 
assumed to have the same total and static conditions of 
temperature and pressure as the primary flow: 


Since ai is uniquely defined by equation (B24), it 
cannot be changed to maintain the maximum clearance 
weight flow at values of (^'/^'cr)|,i greater than 
(K/ Kcr)j Clearance weight flow therefore decreases 
above (K/ Kcr), , . 

From here on, all conditions at station 1 that are the 
same for both the primary and clearance flows (7, 7', p, 
7'i K and p) will have the subscripts p or cl deleted 
and will be identified only by the subscript 1. 

Station 0, - Conditions at station 0 are evaluated from 
the specified geometry, inlet angle, and gas total 
conditions and by using mass continuity between stations 
0 and 1 . The ratio ( V/ is determined iteratively from 


2rj 

h. 

II 

(B14) 

,cl ~ ^1 ,p 

(B15) 

II 

(B16) 

7i,cl ~7i,p 

(B17) 

This implies that 


^l,cl,cr ~ ^Up.cr 

(B18) 


(B19) 

/ V \ / V \ 


% 

II 

(B20) 

/ V \ f ^ \ 


V ^cr / 1 cl,mx ^ l,p,mx 

(B21) 

Pl,d ~ ^l,p 

(B22) 


The flow angle sirid the flow rate wi are 

determined iteratively from conservation of moment of 
tangential momentum between stations 0 and 1. The 
tangential component of velocity is obtained from 


w,l,cl 


''1 


(B23) 


(Wi p+ Wi ci)/?7q 


1 

1— 

1 

O 


1/(7- 1) 

Aq cos ^ j \ 


(B26) 

Also 





(B27) 



(B28) 

Po 

RTq 


(B29) 

( v \ (lyRg.nY^ 


(B30) 

^0 sin ao 


(B3I) 

K,0= ^0 cos ao 


(B32) 


Station 2. -The total temperature at station 2 is 
assumed to equal that at station 1 


7i=7f 


(B33) 


The clearance flow angle and flow rate thus are 


»| cl = sin “ ^ 



(B24) 


CLFRl cos ai^cl (B25) 


and tangential momentum is conserved between stations 
1 and 2. The tangential velocity is a mass-averaged value 
of the primary and clearance flows at station 1 : 

y, 

y 2= r (^^ 1,0 sin ai,„+w-i ,.| sin ai,d) (B34) 
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For no clearance flow this reduces to 

y^^2=^i sina,,p=K^j (B35) 

By using an initial guess for 2> the conditions at station 
2 are evaluated iteratively as follows: 




(B36) 

{ 

yi 

(B37) 

\ycr/2 

^2yRg,Tiy/2 

7’2=7^2|^1 

7-1 / K y] 
7+1 Vl'crAj 

(B38) 

P2=P{- 


(B39) 

P2=P2 1^1 

■ 7+1 Ucr^J 

(B40) 

^2 

RT2 


(B41) 

«2 == COS - ^ 

(W 

(B42) 

W2 = P2^2^2 cos «2 

(B43) 


The assumed K;. 2 is updated until continuity is satisfied; 
that is, 


(B44) 

Station 3, - The total pressure drop between stations 2 
and 3 is obtained from the laminar and turbulent equa- 
tion for flow between parallel plates by using an average 
flow angle between stations 2 and 3. Furthermore, 
moment of tangential momentum is conserved 
(^2^u,2 = From an initial guess of a 3 an average 

turning angle, a hydraulic diameter, a Reynolds number, 
and a flow path length are calculated: 


a2_3 


«2 + «3 
2 


(B45) 


D^ = 2h] cos «2-3 


(B46) 


P 

(B47) 

, _ 

^ ^ cos «2-3 

(B48) 


Friction factor /is defined as follows: 

/=0.24 forRe2<100 (B49) 

f=P~ for 100 < Rc 2 < 3000 (B50) 

/=3.3368xlO-7 Re^-2596 for 3000 <Rc 2 < 3700 (B51) 

/=0.0014+^^ for3700<Rc2 (B52) 

Re2' 

The change in total pressure and the total pressure at 
station 3 are obtained from 



(B53) 

11 

1 

(B54) 

Assuming that 


II 

(B55) 

'*2 = ^1 

(B56) 

and 


V ,= ^2/'-3 
Vua 

(B57) 


all other conditions at station 3 are evaluated from 
P3=Pi(^]) (B60) 
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II 


(B61) 

,/ ^r.2P2'^2 


(B62) 

Finally, a new exit angle 



03-,a„ .(^) 


(B63) 

is calculated and compared with the previously used 
value. This procedure is repeated until successive values 
of are within a specified tolerance. 

Relative conditions at station 3 are calculated from the 
following equations: 



(B64) 



(B65) 



(B66) 



(B67) 

P3=P3(^j 


(B68) 



(B69) 



(B70) 

The rotor incidence angle (from optimum) is calculated 
as follows: the optimum 2 is obtained from the 

Wiesner slip factor correlation (ref. 4), which, when 
adapted to a turbine geometry, takes the form 


where 


^4, hub ■*" ^4, tip 
M,av “ 2 

(B72) 

1 

p(8.16cos ofb|/A/bij 

(B73) 

The last term in the numerator of equation (B71) is a 
correction factor that is applied only if the ratio ^*4 ay /^3 
greater than ejinv Then 

^w,3,opt “ ^w,3,opl ^3 

(B74) 


(B75) 

03 ~ 03, opt 

(B76) 

The loss in total pressure due to incidence is calculated 
as follows; the incidence loss is modeled as in reference 2 
by 

W^{ \ - cos^a^ i) 

2g^ 

(B77) 

where the exponent n is given by the input parameters 
PNEGI and PPOSI for negative and positive incidence, 
respectively. The loss due to incidence can also be 
expressed as a difference between ideal and actual 
velocity: 

Wj~ -cos'^a3j 

) (B78) 

Solving for W'^ . gives 


■= C0S'^Of3 

(B79) 


The relationship 


^1^,3, opt ~ ^3 


( 


1 -Vcos abi/(Nbi)0-7][l - [ (/•4,ay/^3-^lim)/(l -t|im)]^] 


1 - tan a^j/tan 


) 


(B71) 
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^2 (P^yy-"'-’ 

3,1 Y-1 |_ \Pii/ 

implies that 

Pii 

r (7-i)>r2 

*“ lyg^ri 


(B80) 


(B81) 


Substituting equation (B79) for . then gives 


station 4 P 4 is obtained from /?4 by considering the 
combined effects of incidence and rotor viscous loss. 

The value of maximum flow rate per unit area 
(p ^4,mx* well as the value of ( W/ where it occurs 
calculated by subroutine MAXVAL. 
The total relative pressure at station 4 is formulated as 
follows: the incidence loss is modeled by (see previous 
section) 




B^2(i -cos'^a3^/) 


(B86) 


— (B82) 

lygcRT^ 


Station 4. - Flow conditions at station 4 are calculated 
by considering the effects of radius change, incidence, 
and viscous loss. By assuming an initial value of static 
pressure at the hub sector the pressure at each successive 
sector is obtained from simple radial equilibrium. The 
total mass flow of all sectors is then compared with the 
mass flow at station 1, and the next guess of hub-sector 
static pressure is changed accordingly. This is repeated 
until mass flow has converged within a given tolerance. 
The procedure is given here (with calculations performed 
for each sector at the midpoint radius). 

The changes in total temperature, critical velocity, and 
ideal relative total pressure due to radius change are 
calculated from 



(B83) 

W^cr,4=f^cr,3(^J (B84) 

P4,.d=P3(^.j (B85) 

Note that the ideal relative total pressure at station 4 is 
based on the relative total pressure at station 3 and not on 
the relative total pressure due to incidence loss At 

station 3, was calculated for information purposes 
only. As will be shown, the relative total pressure at 


and the loss due to viscous friction in the rotor (rotor 
blade profile, endwall, and secondary flow losses) is 
modeled as in reference 1 by 




k { w \ cos2a3_y+ 


(B87) 


where K is the specified rotor viscous loss coefficient 
(input parameter XKR). The sum of these two losses is 
expressed as the difference between the squares of an 
ideal and an actual velocity to give 


^ 4 = ^^ 3(1 -COS"a3_,) 

+ A'(W"2 cos2«3 ,+ (B88) 


This equation can be solved for W 4 



w^{\ -cos^q:3 -KWl: i 


\^K 




1/2 


(B89) 


with obtained from the relationship 


-A 

4, id 


W^cr.4^'^''' j / P4 

\Pl\d) 


(B90) 


For a given value of static pressure the relative velocity is 
thus calculated by equation (B89). Also, 


T4=T^~ 


(y-iWj 

2yPgc 


(B91) 
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. _ P4 
RTi 

„ /riXT/o-i) 

P4=P.[jr) 


(B92) allowed to exceed 1.0. If it does, the value of 

hub-sector static pressure is adjusted until 4 equals 
1.0. That value of hub-sector static pressure is also 
printed out. 

(B93) jhe absolute conditions are calculated from the 
relative conditions by 


The static pressure at each radial sector above the hub 
sector is obtained from 


</p_ P4(A:-l)[[K(Ar-IV4+ y{k)^_4]/2]^ 
dr g(-r(A:-l )4 


P4(^) =PA^k- 1)+ ^Ar 


(B95) 


where k designates the given sector and (/:-!) designates 
the sector immediately below it. 

For less than the flow follows 

the specified exit angle 13^; for (iv7iv^r\ greater than 
the exit angle is decreased to maintain the 
maximum flow rate. The reduced exit angle is calculated 
from 


04 = COS - 


( P ^) 4 .mx 

_P4W^4^5(1 -BL4) 


(B96) 


where BL4 is the ratio of the blocked area at station 4 to 
the free flow area at station 5. Weight flow and velocity 
components are calculated from 


^4 = P4^4^si^ - BL4) cos 04 

(B97) 

y^uA~ ^4 1^4 

(B98) 

^x,4 = B ^4 cos 04 

(B99) 


If the rotor is fully choked (all sectors), the program 
solves for that value of stator {V/ where rotor choke 
occurred. Also, the maximum value of station 4 hub- 
sector static pressure for rotor choke is calculated and 
printed out. Once the rotor has fully choked, the hub- 
sector static pressure is incrementally decreased by a 
specified amount (input parameter P4RDPC) until limit 
loading is achieved (axial Mach number equal to 1.0). 
The axial component of Mach number is calculated from 


^1/2 



cos 04 


(BlOO) 


yu,4=^uA+^4 

(BlOl) 

II 

(B 102 ) 


(B103) 


(B104) 


(B105) 

P4 

(B106) 

Station 5. - At station 5 the flow area is divided into 
the same sectors as at station 4. All calculations shown 
below are performed for each sector at the midpoint 
radius. The total temperature and tangential momentum 
at station 5 are assumed to be equal to those at station 4: 

II 

(BI07) 

II 

(B108) 

Equation (B107) implies that 


II 

(B109) 

From an initial guess of a 5 the flow 
established iteratively as follows: 

conditions are 

Sin «5 

(BllO) 

Vf^crA ycr,5 

(Bill) 


(B112) 

, , /, K5C0S«5\2/®-^^4»^^\ 

(B113) 
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(B124) 


P5=P5 


‘[-(sj)K):] 


21 7/(7 -1) 


(B114) 


TRQ = 


Wtot>/ A/t6.5(60) 
27rfl 


n - £L 

RT^ 


Kx",4P4^5(^ •" BL4) 


(B115) 


(B116) 


The change in total enthalpy across the rotor (station 3 to 
station 5) for each sector is calculated from 


A/i3_5 = 

and its 


mass-averaged value, from 


(B125) 


A new exit angle is then calculated and compared with the 
previous value. 

«5 = tan-i(^^^ (B117) 

The procedure is repeated until successive values of are 
within a specified tolerance. 

Performance parameters, -Tht performance param- 
eters are calculated as follows: for each sector the static 
and total ideal enthalpy change from station 0 to station 5 
is calculated from 

Mass-averaged values of ideal enthalpy change, power, 
and torque are established from 


>^tot= Yi (B120) 

k=\ 

N 

Ah(k)Q_^w(k) 

= — (B121) 

w-tot 


N 

£ Ah' (k)i^^w{k) 

Ahi_^=— (B126) 

Wtot 


The aerodynamic efficiencies for each sector are 
calculated from 


_ A^3-5 

^aed “ Af, 

^«0-5 


(B127) 


^aed ~ 


^^3-5 

A/j6-5 


(B128) 


Disk friction and exit tip clearance loss per unit mass flow 
are obtained from (ref. 3) 


0.02l25p3D|r2 



Tci ~ 

\ — 

^4, tip 


(B129) 


(B130) 


Actual total enthalpy, actual power, and actual torque 
across the rotor are 


N 


Ao-5 = ■ 


Yj Ah' (k)Q^^w(k) 
k=\ 


VV| 


tot 


^^6-5 

550 



^^act ~ ^3-5 ^df ^cl 

(B131) 

(B122) 

p '^tot^^ ^^act 

“ 550 

(B132) 

(B123) 

T~D/^ ^^act 

2rrfi 

(B133) 
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and actual overall efficiencies are 


/?std = 287.05 J/(kg-K) (53.35 ft-lbf/lbm-°R) 


(B146) 


r,= 

A/»o-5 

(B134) 

r,’= 

^^6-5 

(B135) 

Mass-averaged static and total 
calculated by 

temperatures are 

X; T{k)^w{k) 

•^to. 

(B136) 

,\ 

T'{k)^w{k) 

^tot 

(B137) 

Static pressures at the hub and shroud and average static 
and total pressures are obtained from 

/dp\ _p(1)5K2(1)„5 

\ dr /hub «f^(l)4 

(B138) 

/dp\ _p(N)5|/2(N)„5 

\^^''Ahr gc-r{N)4 

(B139) 


(B140) 

(ILt 

(B141) 

„ P5,hub “*“/^5,shr 

P5,a.~ 2 

(B142) 

/fs\-ydy-\) 
P5.av=PS,av\^Y J 

(B143) 


Equivalent parameters are calculated on the basis of the 
following standard conditions; 

rsid = 288.2 K (51 8.7 °R) (8144) 

/?s'td = 101 325 N/m2 (21 16.22 lbf/ft2) (8145) 


7std “ 1 

(B147) 

The equivalent parameters (flow, speed, work, power, 
torque, total-to*static pressure ratio, and total-to-total 
pressure ratio) are 

^equ - ^ 

(B148) 

N - -AL 

(B149) 

WRK = 

yy ivrvgqu ^ 

(B150) 

p _ ^aci^ 

^equ - 

6V6^cr 

(B151) 

TRO - 
^ ^Vequ “ ^ 

(B152) 

/ Po \ — F] (7std~0*^ ^^0-5 

\/^5/equ L ^crTstd^std^std 

-7sid/'(7sKl - 1) 


(B153) 

- fl (fstd- IV AAo_ 5 
\/^5/equ L ^crTstd^std ^std 

~ 7sid^(7s[d ~ t ) 

(B154) 


where 



(B155) 

8=^ 

Pkd 

(B156) 

/ 2 \ 7std'^(7std ~ n 

''*'‘'V7s.d+l/ 

/ 2 y/(7-i) 
^\7+ 1 / 

(B157) 
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The overall dimensionless parameters of specific speed, 
specific diameter, work factor, and loading coefficient 
(blade-jet speed ratio) are defined as follows: 



(B165) 


^5,av 


P5,av 

^^5,av 


(BI58, e=(;-^|)(T^)' 


(B166) 





(g^ A/z6-5)«-75 


(B159) 




(B160) 




A/?i_5 




(B161) 


^3 

{Ig^ Aho-sV^^ 


(BI62) 


Subroutine MAXVAL 

Subroutine MAXVAL finds the maximum value of a 
function over an interval. The chosen function must be of 
a form such that a maximum exists within the interval 
(other than at the end points). MAXVAL is used to 
calculate the maximum values of flow rate per unit area 
at the stator exit (pL) i and rotor exit (PHO4 sectors. The 
actual mass flow calculations are performed by the called 
subroutines ROVST and ROVROT for the stator and 
rotor, respectively. Given the interval, the function to be 
maximized, and an initial number of divisions of the 
interval, MAXVAL determines the function’s maximum 
value and its location in the interval. 


Subroutine ROVST 

Subroutine ROVST calculates the value of density 
times velocity at the stator exit (station 1) for the assumed 
stator loss model and a specified (K/K^.j.)j. The formula- 
tion is as follows: 




^ 2yRg,T{ y'2 


(BI63) 



(B164) 


If the total pressure at station 1 is calculated by specifying 
IOPTS= 1 (specifying the pressure ratio p{/po). 



(B167) 


If the total pressure at station 1 is calculated by specifying 
lOPTS = 2 (specifying ^30), 


Pi =Po 


4 


1 - g3D - Q 
(l-e3D)(l-0 


] 


7 At “ I ) 


(B168) 


Finally, 


P\ 




7/(7- 1) 


(B169) 


Pi 


P\ 

RTx 


(B170) 


pK-piLj 


(B171) 


Subroutine ROVROT 

Subroutine ROVROT calculates the value of density 
times velocity at the rotor exit sectors (station 4) for the 
assumed rotor loss model and a specified ( W/ B^cr)4* The 
equations used are as follows: 



7'4=T4'- 


^ygc^ 


(B172) 

(B173) 


The static pressure is calculated from equation (B88), 
which expresses the sum of the two rotor losses (incidence 
and viscous friction) as the difference between the 
squares of an ideal and an actual velocity 
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Then 


+ A:(fr2 cos2a3,i+ (B174) 

By using 




(B176) 


equation (B174) can be rewritten as 


Pa 




7/(7 - 1) 


PA = 


Pa 


RT4 
PV=PA^A 


7/(7- i) 


Pa ~ P4,iti 



/ IV^ 

V 

COs2„, j - ( ^^2 J (1 -COS"0,,.,) 

( 

'7+l\ 

.7-1/ 

_ / 1^4 y 


(B178) 

(B179) 

(B180) 


(B177) 
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